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Abstract 

We re-examine the justification for the imposition of regular boundary 
conditions on the wavefunction at the Coulomb singularity in the treat- 
ment of the hydrogen atom in non-relativistic quantum mechanics. We 
show that the issue of the correct boundary conditions is not independent 
of the physical structure of the proton. Under the physically reasonable as- 
sumption that the finite size and structure of the proton can be represented 
as a positive correction to the Coulomb potential, we give a justification for 
the regular boundary condition, which, in contrast to the usual treatments, 
is physically motivated and mathematically rigorous. We also describe how 
irregular boundary conditions can be used to model non-positive correc- 
tions to the Coulomb potential. 
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1 Introduction 



One of the unsatisfactory features of Old Quantum Theory was the primacy of 
the concept of 'quantum number', which does not arise from more basic physical 
principles thereby leaving a certain mystery as to the true origin of discreteness 
in quantum theory and the distinguished nature of the positive integers. It was 
to remove this limitation that Schrodinger embarked upon a programme of re- 
casting quantisation as an eigenvalue problem, studying the hydrogen atom as 
his first example [|I|. 

Of course, eigenvalue problems such as those encountered in quantum theory 
require the specification of boundary conditions in order to be well-posed. In the 
hydrogen atom case, there is a particular issue about what boundary conditions 
should be imposed on the wavefunction at the Coulomb singularity. This was 
understood by Schrodinger, who imposed that the wavefunction be finite. With 
this condition, the negative energy eigenvalues of the idealised Coulomb Hamil- 
tonian precisely replicate the Bohr levels, a success which confirmed the validity 
of the new approach. As Schrodinger wrote 

The essential thing seems to me to be, that the postulation of "whole 
numbers" no longer enters into the quantum rules mysteriously, but 
that we have traced the matter a step further back, and found the 
"integralness" to have its origin in the finiteness and single- valuedness 
of [the wavefunction] 

But is finiteness of the wavefunction an axiom of quantum mechanics? Al- 
though this issue is often considered in quantum mechanics textbooks, there does 
not appear to be a satisfactory treatment in the literature - see for a discus- 
sion of the usual treatments. In the present paper, we examine this problem in 
detail for the particular case of the hydrogen atom. One of our main aims will be 
to provide a physically motivated and mathematically rigorous justification for 
the condition that the wavefunction should be finite at the Coulomb singularity. 
We shall also consider the interpretation of other possible boundary conditions as 
models for nuclear structures differing greatly from that of the physical proton. 

Our general methodology, which is discussed at length in p[, starts from the 
distinction introduced in between true problems and idealised problems. In a 
true problem, interactions are modelled by a smooth potential in the Schrodinger 
equation, reflecting our expectation that only smooth potentials occur in nature. 
In the case of the hydrogen atom, one can identify a class of true problems 
consisting of models in which the finite size and structure of the proton is modelled 
by a smooth potential which deviates from the Coulomb form within the nuclear 
radius, 'rounding off' the Coulomb singularity. For simplicity, we shall mainly 
consider spherically symmetric potentials of this form, but we make no other 
restriction on our class of true problems. For true problems, we will take finiteness 
and smoothness of the wavefunction to be axiomatic for quantum mechanics. 
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Idealised problems represent limiting cases of true problems. As such, they can 
possess distinguished points at which some of the structure of a more fundamental 
true problem has been simplified. In the case of the hydrogen atom, the Coulomb 
singularity is a distinguished point, for it represents the physical proton of the 
true problem. Note that the idealised problem need not be singular at such 
points. For example, a true problem in which the potential changes steeply but 
smoothly from one value to another could be idealised by a potential step. One 
might idealise a true problem with a potential compactly supported within a 
small neighbourhood of the origin, by an idealised problem with no potential 
at all, but with the origin as a distinguished point (see [0], where this case is 
treated in depth). At such distinguished points in idealised problems, boundary 
conditions can only be physically justified by reference to the true problem one 
intends the idealisation to represent. 

Accordingly, in Section 2, we determine the class of allowed idealised problems 
for the case in hand by considering limits (in the strong resolvent sense ||15|| ) 
of sequences of Hamiltonians for true hydrogen atom problems whose 'nuclear 
radii' tend to zero. We thereby determine a 1-parameter family of idealised 
problems with different boundary conditions at the Coulomb singularity, all of 
which correspond to the limiting behaviour of true problems from our class. The 
spectra of these idealised problems can differ markedly in the S'-wave, although 
for higher angular momenta, all idealisations yield the Bohr levels. 

Clearly, in order to justify the choice of regular boundary condition at the 
origin (the unique choice leading to the Bohr levels in the S'-wave) we must 
make further restrictions on our class of true problems. In Section 3, we make 
the physically reasonable restriction to true problems whose potentials represent 
positive corrections to the Coulomb potential localised within a small nuclear 
radius. Under these assumptions, we prove that the true problem has energy 
levels which differ only slightly from those of the idealised system with the regular 
boundary condition (i.e. the Bohr levels). There are, of course, well known 
estimates from first order perturbation theory of the finite size corrections to the 
Bohr levels under the further assumption that the nucleus is a smooth positive 
charge distribution. However, our treatment will cover a wide class of nuclear 
models in a non-perturbative fashion. We also note that this argument is not 
usually given as the 'correct' understanding of the point in question. We will also 
show that there exist true problems with arbitrarily small nuclear radius (with 
non-positive corrections to the Coulomb potential) whose S'-wave energy levels 
differ greatly from the Bohr levels. 

We also provide a second justification by considering limits in the sense of 
strong resolvent convergence of sequences of true problem Hamiltonians whose 
nuclear radii tend to zero, intepreting this as a means of 'regularising' the point 
charge. By definition, in the general case, the class of possible self-adjoint limits 
is precisely the class of idealised problems, which (as we mentioned above) is a 
1-parameter family of operators with different boundary conditions at the ori- 
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gin. However, restricting to regularisations which can be modelled by positive, 
compactly supported corrections to the Coulomb potential, we will prove that 
all such regularisation schemes converge rigorously to the idealised Hamiltonian 
with regular boundary conditions. 

In Section 4, we turn to the intepretation of irregular boundary conditions 
at the Coulomb singularity, by developing a scattering length formalism for sys- 
tematically matching any given true problem to its 'best fit' idealisation. This 
formalism also plays a role in the proof of our convergence theorems mentioned 
above. In Section 5, we examine the situation for angular momenta i > 1. We 
find that the Bohr spectrum in these sectors is more stable against finite size 
effects than in the S'-wave. Section 6 contains the rigorous proofs of some of the 
results stated in earlier sections. 

One of the motivations for this work was the recent study undertaken by 
Kay and the author of the role of model dependence in systems which interact 
with small objects ||^, and in particular the 'principle of sensitivity' introduced 
in 1^ (see also Q]) which relates the range of possible large scale (low energy) 
behaviour of a class of true problems to the range of possible behaviour exhibited 
by the corresponding class of idealised problems. We conclude in Section 7, by 
discussing the relation of our present results to the principle of sensitivity, with 
which we find good agreement. 

2 True and Idealised Problems 

The Hamiltonian corresponding to the idealised hydrogen problem is given as a 
differential operator by 

//ideal = -A + ^, (2.1) 

r 

where we employ units in with ft = Aneo = 1 and in which the reduced mass of 
the electron is set to | (7 is negative). The Bohr radius is thus given by ao = 
2/I7I. Our Hilbert space of wavefunctions is 7i = L^(R^,c?^r) with inner product 
denoted (■ | ■). To define i/ideai rigorously on H, we must specify its domain, 
which inevitably begs the question of the appropriate boundary conditions at 
the Coulomb singularity. In order to circumvent this, we consider the class of 
true problems, in which the singularity is rounded off within some nuclear radius, 
and for which we assume finiteness of the wavefunction as an axiom. We then 
define the class of idealised Hamiltonians to consist of the self-adjoint limits in the 
strong resolvent sense of sequences of true problem Hamiltonians whose nuclear 
radii tend to zero. 

Let V be the class of all measurable real-valued functions V{r) on R"*" com- 
pactly supported within some radius of the origin such that r + ^(|r|) 
is smooth on R^. Then the class of true problems consists of all Hamiltonians 
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defined as the closure of an operator of form 

- A + - + V(r) on C^(R^) C H (2.2) 
r 

for V ^ V. We refer to the radius of support of V as the nuclear radius. Opera- 
tors of form ( p.2|) are automatically essentially self-adjoint as a consequence of the 
Kato-Rellich theorem by hypothesis on V (see Theorem X.15 in JTB]). This defi- 



nition encodes two desirable features of true problems: firstly that only smooth 
potentials occur in nature, and secondly that the wavefunction is everywhere fi- 
nite. The latter follows from our choice of domain, which has already injected 
regular boundary conditions at the origin. 

We now define the class S of all sequences of true problem Hamiltonians 
Hn = —A + 'j/r + Vn{r) whose nuclear radii (i.e. the radii of the supports of the 
Vn about the origin) tend to zero. Our idealisations are the self-adjoint limits of 
such sequences in the strong resolvent sense. 

The following result, which we prove in Section 6, classifies the idealised 
Hamiltonians as the self-adjoint extensions of the symmetric operator -ffideai on 
the domain C~(m3\{o}). 

Theorem 2.1 Let //ideal = -A + 7/r on Co~(M^\{0}). (a) Let {H^} E S. // 
{Hn} has a self-adjoint limit H in the strong resolvent sense, then H is a self- 
adjoint extension o/i/jdeai- (b) Furthermore, all self-adjoint extensions of H-^^^^x 
arise as the strong resolvent limits of sequences in S. 

We note that these extensions have also been exhibited as limits of sequences 
of self-adjoint operators by Albeverio and co-workers ||^, where however, the 
convergence is in the norm resolvent sense. We motivate our choice of convergence 
by two observations. Firstly, the strong resolvent convergence of a sequence of 
self-adjoint Hamiltonians to a self-adjoint limit is equivalent (Trotter's Theorem 
[|T5[] ) to the statement that 

giH„t^ ^ (,iHt^ ^jj ^ g ^_ (^2.3) 

Thus strong resolvent convergence is a natural candidate for the notion of dynam- 
ical convergence. Under certain circumstances, it can be shown that if if„ — H 
in the strong resolvent sense, then the corresponding M0ller wave operators 
VL^i^Hn, Ho) converge strongly to fl'^^H, Ho), where Ho is a suitable comparison 
dynamics f^. 

Secondly, the spectrum can contract in the limit of strong resolvent conver- 
gence (although it cannot expand) in contrast to norm resolvent convergence, in 
which the spectrum can neither contract nor expand. (See Theorem VHI. 24(a) 



and the following discussion in |]T5|-) This property of strong resolvent conver- 
gence will be important in the sequel; furthermore norm resolvent convergence is 
too strong for our purposes: one can see that there are no analogues of some of 
our results in this sense of convergence. 
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We also note that the results in concern sequences of scaled Hamiltoni- 
ans, whereas (due to our weaker notion of convergence) we are not restricted 
in this way. In addition, the scattering length formalism developed in Section 
4 provides a transparent physical interpretation for our convergence results. In 
these differences are discussed at greater length in the case of short range 
(non-Coulombic) potentials. 

We quickly review the classification of the self-adjoint extensions of -ffideai 
along the lines of the discussion in . The deficiency indices [0 of -f^idcai are 
defined to be the dimensions of the deficiency subspaces = keT{H*^^^i =f i), 
i.e. the number of independent solutions to 

i/rdeal^ = ±#- (2.4) 

ifj^gg^i may be seen to act 'classically' (cf. Proposition 2 in the appendix to section 
X.l of [|16]) i.e. such that any such solution solves 

(^-A + ^^tp = ±iip (2.5) 

as a differential equation, and so usual techniques of partial differential equations 
suffice. In the familar way, we use the decomposition 

oo 

L2(m^ dh) = L2(m+, r^dr) (g) ICe (2.6) 

£=0 

where Ke is the subspace of L^(S^,(if2) spanned by Y£,_£, . . . , Y^,^, to separate 
over the basis of spherical harmonics. We make the unitary transformation 
U : L^(K"'", r^rfr) L^(M"'", rfr) given by {Ux){f') = ''"x(^) to obtain the radial 
Hamiltonians for each angular momentum sector 

which we define on C^(0, oo) c L^(E"'", dr). (Studying the hi on this domain is 
actually equivalent to studying if ideal on V, the set of finite linear combinations of 
terms of form /(r)Y^^m(^, V^) where /(r) G C^(0, oo). However, this can be seen 
to yield the same deficiency subspaces and self-adjoint extensions as are obtained 
with V = C^{R^\{0}).) 

Solving h}uf = ±iuf as an ODE yields, as solutions square integrable at 
infinity 

uf = W_,^/2(±,)i/2,,+i/2(-2^(±^)^/V) (2.8) 



where Wk-Jz) is a Whittaker function 18 . For i > 1, these functions are not 



square integrable at the origin and so the hi are essentially self-adjoint. In the case 
£ = 0, however, the functions Uq are square integrable, so ho has deficiency indices 
(1,1) and a 1-parameter family of self-adjoint extensions labelled by f/(l), or 
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equivalently (under an obvious correspondence) by the extended real hne ]RU{oo}. 
These may be denoted[] where L G M U {00} and has domain 

D{h^) = {ue L2(m+) I m, m' G ACiociR^y, u" + ir-^u G L2(m+), 

[L-i + 7(^(1) + ^(2))]no + Ml = 0} (2.9) 

where uq = hmr^o+ '"('^)! ""i = hm.r^o+ ''^~"'^('"(^) ^ ""0(1 + T'^logdTk)))- Here 
AC\oc{^) denotes the set of locally absolutely continuous functions on the 
positive half-line, and '^{z) = d/dz\ogT{z) is the di-gamma function HTs). Ui 
can be regarded as a 'Coulomb modified' first derivative (from the right) of u{r) 
at r = 0. Note that the usual Hamiltonian, with regular boundary conditions, 
arises from the choice L = 0. 

The spectrum of Hq in the case 7 < was first discussed by Rellich (see 
for the 7 > case). For each value of L, there are infinitely many non-degenerate 
negative eigenvalues of h^. The eigenfunctions are precisely those solutions v{r) 
to 

K^v (2.10) 




which are square integrable and lie in the domain (|2.9| ) above. The first require- 
ment entails that v{r) is given (up to normalisation) in terms of the Whittaker 
function by 

^W=Wm/2.;|(2«^^)- (2-11) 

Defining vq = limr^o+ v{r), Vi = limr^o+ r~^{v{r) — vq{1 + •yr log(|7|r))), we find 
that 

— = j log ^ - ^(1) - ^(2) + ^(1 + 7/2k) - k/j (2.12) 

^0 L I7I 

and hence (by our second requirement) that eigenvalues occur for E = — 
satisfying 

where the function G is defined by 

G{z) = ^(1 + z)- log 1^1 - l/{2z). (2.14) 

G{z) is everywhere increasing in 2; G M except at its poles, located aX z = —n for 
n = 0, 1, 2, . . ., and so G{z) increases from —00 to 00 between each consecutive 
pair of poles (see Fig. 2.1). There is therefore precisely one zero Zn of G{z) be- 
tween each consecutive pair of negative integers Zn G (— n, —n + 1) for n = 1, 2, . . . 
and precisely one solution of ( |2.13| ) for 7/ (2(— £^)^/^) in the interval (zn+i, ^n] for 
each n = 1, 2, . . .. We label the corresponding energy by En- For L > 0, there is 

^Our parametrisation differs from that in H. This wiU be convenient later. 
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also an eigenvalue Eq, with 7/(2(— i?o)^''^ ^ i^i^O). In addition, the special case 
L = - which corresponds by ( p.9|) to the case of regular boundary conditions 
at the origin - gives rise to the usual Bohr levels 

En = -^, n=l,2,.... (2.15) 

However, we see that by changing L, we can radically alter the eigenvalues. The 
eigenvalues En for n > 1 are subject to maximum fractional changes of order 
1/n, and for all L > 0, there is an additional eigenvalue £"0 lying well below the 
usual Bohr levels. Thus the ground state and first few excited energy levels can 
be quite sensitive to the value of L. However, if \L\ -C ao, we can see (comparing 
with Fig. 2.1) that the energy levels are little changed from the Bohr levels. 
Indeed, in such cases, because \1/(1 + ^) = —{C, + n)^^ + 0(1) as ,^ ^ —n — 1 for 
all n = 1, 2, . . ., we may write the second equality in ( |2.13| ) in the approximate 
form 

(2.16) 



1 1 



7L 7/2(-E)i/2 +rz 
obtaining the approximate solutions 



7^ 7^ 



Thus 2L/ao plays the role of the Rydberg correction (or quantum defect) of 
atomic physics ||T3|. In addition for < L <C ao, because G(^) = — (2^)~^ + 
0(log 1^1) as ^ ^ 0, we have an energy level Eq lying well below the Bohr spectrum 

Eo ^ (2.18) 

The presence of such a deeply bound state is initially puzzling, and we shall 
return to this in Section 4. 

As noted above, for i > 1, the operators hg are essentially self-adjoint on 
the domain C^(0, 00) and so there is a unique self-adjoint extension hi whose 
eigenvalues are simply the familiar Bohr levels: 

^'^'■--M^^ «=1,2,.... (2.19) 

By analogy, we will abuse notation and write Jiq to denote Hq. We may therefore 
assemble the radial Hamiltonians to give the full Hamiltonian (recalling the 
decomposition (|2.6|)) defined by 

00 

= U*h^U ®l®^U*heU (^1 (2.20) 
e=i 
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which acts on the S'-wave as /iq and as hi on the £'th sector for ^ > 1. The 
spectral properties of are immediately given by the discussion above. 

The case corresponds to the usual choice of idealisation with regular 
boundary conditions. This operator is the closure of if ideal on C^(M^) and would 
therefore result if we treated the idealised problem in the same way as true prob- 
lems by assuming finiteness of the wavefunction as an axiom. We emphasise that 
Theorem shows that it would be inconsistent to do this. 

To summarise, there is a 1-parameter family of idealised Hamiltonians arising 
as self-adjoint strong resolvent limits of sequences of true problem Hamiltonians 
whose nuclear radii shrink to zero. The idealisations are the self-adjoint ex- 
tensions of — A+7/r on C^(m^\{0}) which are labelled by a single real parameter 
L G MU {c>o} (later, we will identify L as a Coulomb modified scattering length). 
The idealisations differ from the usual Hamiltonian H*^ only in the S'-wave, in 
which the low-lying energy levels can exhibit significant fractional changes from 
the Bohr levels, depending on the value of L. In addition, for all non-zero values 
of L, there is an additional bound state lying below the usual Bohr spectrum. 
We will discuss this bound state later in Section 4. 

3 The Regular Boundary Condition 

We now turn to our rigorous justification of the regular boundary condition. 
Central to this is the mild physical requirement that the correction potential V is 
positive and has nuclear radius a <^ a^. Under these assumptions, we will see that 
the spectrum of the true problem is well approximated by the Bohr levels. In the 
following, for any self-adjoint operator H bounded from below, fin{H) denotes 
(when it exists), the n'th discrete eigenvalue of H, counting with multiplicity and 
in increasing order. 

Given positive G V, we define the S'-wave radial true Hamiltonian /io,true = 
ho + V. Suppose V has nuclear radius a. Because V is positive, the n'th eigen- 
value /in(/io,true) IS bouuded below by the n'th Bohr level fin(ho). Moreover, it 
is resonable to suppose that an upper bound should be provided by the n'th 
eigenvalue of /iq, the Hamiltonian describing an impenetrable nucleus of radius a 
which corresponds to the heuristic potential 



In fact, this turns out to be the case. We define /iq defined on L'^{[a, 00), dr) by 




(3.1) 




(3.2) 



where 



V' = {i;\{r- a)-V e C^i[a, 00))} C L^{[a, 00), dr). 



(3.3) 
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The domain V is chosen to ensure that all eigenf unctions vanish at the radius 
r = a. In the Appendix, we prove that this Hamiltonian has precisely one 
eigenvalue between each consecutive pair of Bohr levels and no other eigenvalues. 
(In addition there is essential spectrum [0, oo).) Note that this Hamiltonian is 
defined on a different Hilbert space from /io,truo- In Section 6, we prove 

Theorem 3.1 Let /io,true = + V for some positive V &V with nuclear radius 
a > and define h'^ by ( p'.i^ j for this choice of a. Then Jiq, /io,true and H'q are 
hounded from below and possess infinitely many eigenvalues satisfying 

finijlo) < yU„(/lo,truc) < f^niho) (3.4) 

for each n = 1, 2, . . .. 

In the case a -C oq, the energy levels of H'q are estimated in the Appendix as 



Hence (p.4|) may be re-written 



(3.6) 



and so, using a/oo ~ 10~^, the maximum fractional error for the physical hydro- 
gen atom is of order 10~^. 

Hence we see that the energy levels of the true problem are close to the Bohr 
levels, provided the correction to the Coulomb potential is positive and that the 



nuclear radius is smaller than the Bohr radius. Note that the bound ( |3.6|) is 
considerably weaker than the usual perturbative estimate of the energy shift due 
to the finite size and structure of the nucleus [|1^], which assumes the nucleus is 



a smoothed out region of positive charge and works to first order in perturbation 
theory, obtaining energy shifts of order (a/oo)^. However, the result above is 
non-perturbative and covers a much greater range of nuclear models, including 
many for which first order perturbation theory would break down. In the next 
section, we will describe how a general true problem can be modelled by the 
idealised problem with potentially irregular boundary conditions, and how the 
usual perturbative result may be derived within this formalism. 

As a second justification of the regular boundary condition, we consider 
sequences of Hamiltonians with positively corrected Coulomb potentials. Re- 
call that S denotes the class of sequences of true problem Hamiltonians Hn = 
—A -|- 7/r -|- Vn{r) whose nuclear radii tend to zero. Our space of general reg- 
ularisation schemes is therefore S. In general. Theorem |2.1| - which defines our 
notion of idealisation - shows that class of the possible self-adjoint limits of such 
sequences in the strong resolvent sense is precisely the class of self-adjoint exten- 
sions of -ffideai on C^(R^\{0}). Thus, in general, the limit (when it exists) is not 
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the self-adjoint extension with regular boundary conditions. However, restricting 
to the subclass S+ C S of sequences iJ„ = — A + 7/r + \4(r) in which the poten- 
tials Vn are positive, it turns out that all such sequences converge in the strong 
resolvent sense. 

Theorem 3.2 All sequences in S+ converge in the strong resolvent sense to , 
the self-adjoint extension o/i/idcai with regular boundary conditions. 

Thus any regularisation scheme in E"*" selects the regular boundary condition, 
providing our second justification. 



4 Scattering Length Formalism 



In this section, we consider how a general true problem may be approximated by 
an appropriate choice of idealised problem. The starting point for our discussion 
is the low energy expansion of Coulomb modified scattering theory developed 
by Lambert [|l^. This will enable us to examine the behaviour of both true 
problems and idealisations at energies small in comparison with the nuclear scale, 
i.e. E -C This regime covers the atomic spectrum, whose characteristic scale 
is the Bohr radius, 5 orders of magnitude larger than the atomic radius for the 
physical hydrogen atom. The 'best fit' boundary condition can then be found 
by matching the low energy behaviour of the idealisation to the true problem. 
The procedure below is the direct generalisation of the treatment of the non- 
Coulombic case in |]^. 

In the S'-wave, we therefore examine the regular solution u{r) to 
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+ - + Vir) 



u 



(4.1) 



where V{r) obeys r^F(r) ^ as r ^ and is compactly supported within some 
radius a of the origin. These conditions allow all potentials in our class V. 

The Coulomb modified phase shift 5o(^) is defined by writing the asymptotic 
form of the regular solution to ( |4.1| ) as 



u{r) ~ sin(fcr - (7/(2A;)) log2A;r + cro(fc) + (5o(A;)) 
where cro(/i;) is given by 



V{l + i-i/2k) 
T{l-i-i/2k)' 



(4.2) 



(4.3) 



Defining 



Tvy/k 



-icot5o(fc) + 1 



2k 



Z7 
2k 



7log^, (4.4) 
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we have the (Coulomb modified) low energy expansion (see e.g. |I2| ) 

Z'^{k) = -j + ^r,e + Oik^) (4.5) 

where L and are the Coulomb modified scattering length and effective range re- 
spectively. This is the analogue of the familiar low energy expansion of scattering 



theory for short range forces |11 



It can be shown |^ that the scattering theory of the idealised Hamiltonian 
is described by the S'-wave low energy expansion 

Z'^ik) = -\. (4.6) 

(For angular momenta i > 1, there is, of course, no scattering.) Thus the operator 
has a low energy expansion which is exact at lowest order. We are also able 
to interpret the parameter L as the (Coulomb modified) scattering length of H^. 
It is therefore clear from this and equation ( [4.5| ) that scattering theory of 
describes the leading order scattering theory of any true problem with scattering 
length L at low energies. Moreover, because the low energy expansion may be 
analytically continued to discuss bound states, should also well-approximate 
the bound states of any true problem with scattering length L at sufficiently low 
energies. It is therefore possible to select the 'best fit' idealisation to a given true 
problem simply by matching the scattering lengths. 

In order to compute the scattering length of a given true problem, it suffices 
to consider the regular solution u{r) to (|4.1| ) at zero energy, i.e. k"^ = 0. Lambert 



showed |]T2| that the scattering length is then given by 

, _ 0(q) 1 (u'jr) 
e{a) e{af \u{r) 




(4.7) 



Because this formula allows us to select the best fit idealisation to the true prob- 
lem, we refer to it as the fitting formula. Here, 



I7l 



1/2 

I ^i(2(|7|r)^/^) (4.8) 



and 

e{r) = -7r(|7|r)i/2jy^^2(|7|r)^/2) (4^9) 

(where Ji{z) and Ni{z) are the Bessel and Neumann functions of order 1 |T^) 
are the regular and irregular solutions to ( [4.1| ) in the special case = 0, = 0. 
Note that (f) and 9 have Wronskian (f)'9 — O'cf) = 1. The effective range may also 
be determined using u, and 6 [T^ . 



It is, of course important to know the range of energies over which this approx- 
imation holds good. In B, this is referred to as the question of heliev ability. As 
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a simple necessary condition - derived by requiring that the second term in ( ^4.5| ) 
be small compared with the first - we say that the idealisation is believable at 
energy E only if 

1^1 < (4.10) 

\Lro\ 

where L and are the scattering length and effective range of the true problem. 

As we noted in Section 2, for L > oq, the S-wave energy levels exhibit signifi- 
cant fractional changes from the Bohr levels. From this it is to be expected that 
the low-lying energy levels of an atom differ markedly from the Bohr levels if the 
nucleus has scattering length of the order of, or greater than, the Bohr radius. 
It is therefore important to know whether the scattering length can take such 
values. In fact a simple argument shows that for any L G MU {oo} there exist po- 
tentials V{r) of arbitrarily small support with scattering length L: let r \—>- ip{\r\) 
be any smooth spherically symmetric function on with ilj{r) bounded away 
from zero for r < a and iplr) = r~^((/)(r) — L9{r)) for r> a. Then the potential 
V{r) defined by 

V{r) = ^-^Ailj-- (4.11) 
r 

is supported within r < a and the effective potential 7/r + V{r) is smooth on 
R^. Moreover, by (|4.11|) and the first property, V{r) has scattering length L. We 
have therefore proved 

Theorem 4.1 For any L G KU {00}, and any a > 0, there exists a potential 
V E V with nuclear radius less than a and with (Coulomb modified) scattering 
length L. 

Thus there exist true problems with 'nuclei' of arbitrarily small radius whose 
energy levels at sufficiently small energy (i.e. for large enough principal quantum 
number) are well-approximated by the energy levels of if^, for any L G RU {cxo}. 
One can use 'believability' to determine the range of energies for which this 
approximation is valid. Provided the effective range (and all higher parameters 
in the low energy expansion) is small, the approximation could include the energy 
levels of low principal quantum number. Concrete examples are given in |^ . 

We now restrict to the case of positive corrections G V to the Coulomb 
potential. We take the nuclear radius to be a and define u{r) and 0(r) as above, 
noting that they satisfy 

(t){r)u"{r) - M(r)0"(r) = V {r)<i){r)u{r). (4.12) 

The first zero of 0(r) away from the origin occurs at approximately 1.8ao, and 
that of 9{r) at ai ~ 0.6ao. Restricting to the case a < ai (which therefore includes 
the physical case a ~ 10~^ao), we have > on (0, a) and we may write 

(t>{rf^ { = (t>{r)u'{r) - M(r)0'(r) = f V {r')(l){r')u{r')dr' . (4.13) 
dr \ (bir) I Jo 
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Because r'^V{r) ^ as r ^ and u, both obey regular boundary conditions at 
the origin, we may divide through by 0^ to yield an integrable function, leading 
to the following integral equation for u{r) 

u{r) = (f){r) + (f){r) [ dr'{(t){r'))~^ [ dr"V{r")(l){r")u{r"). (4.14) 



Since 0(r) > on (0,a), we conclude that u{r) is also strictly positive on this 
interval, for by the second equality in ( [4.13|) , u could vanish only if u' is positive. 
Therefore, because u{r) is initially positive, it must remain so on this interval 
{y{r) is sufficiently regular that u{r) remains bounded on (0,a)). Hence we 
obtain the relation 

n!{rl_<m^ r,V{r')ct>{r')u{r')dr' ^ ^ 
u{r) 0(r) (l){r)u{r) 

Inserting this in the fitting formula (|4.7|) and using the fact that and 9 have 
unit Wronskian and are positive on (0,a), we find that the Coulomb modified 
scattering length for such a potential satisfies 

0<.<^. (4.«) 

The trivial case V{r) = yields vanishing scattering length. We have thus proved 



Theorem 4.2 Let V E V be positive and supported within [0,a), where a < ai. 
Then the scattering length of V is bounded by 

0<L<^= L^,.. (4.17) 

Furthermore, if a < ao/8, then a reasonable approximation to the upper bound 
is provided by a, the nuclear radius, so we can write 

< L < a. (4.18) 



We can substitute this inequality into the approximate expression ( |2.17|) in 



order to re-derive equation ( p.6[ ) as a useful consistency check (although the 
derivation in this section truncates the low energy expansion at first order and 
is therefore less rigorous than our previous derivation). In fact, we show in the 
Appendix that the n + I'st eigenvalue of /ig provides an upper bound on the 



n'th eigenvalue of and hence (by Theorem 3A.) on the n'th eigenvalue of /io.true 
for a positive correction to the Coulomb potential. Note that, for < L < a, 
the idealisation also exhibits a energy level Eq{L) ~ —L~^ lying well below 
the Bohr spectrum, as we saw in Section 2. However, we can use believability to 
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show that this is not a behevable feature of the true system, at least in the case 
where the effective range is of the order of the nuclear radius, when the range of 
believability is 

\E\ «^<^ (4^9) 

SO the bound state Eq ^ — falls outside this range. The presence of this 
bound state is therefore a result of the idealisation process which does not in this 
case correspond to any feature of the true problem. This is reassuring, because 
it is, of course, impossible for a positive correction to the Coulomb potential to 
introduce a new energy level below the Bohr levels. This provides a simple ex- 
ample of the manner in which the components of the scattering length formalism 
and believability work together. 



5 Higher Angular Momenta 

The discussion in Sections 3 and 4 demonstrated two features of the true problem 



in the S-wave. Theorem showed that for positive corrections to the Coulomb 
field of small nuclear radii, the true spectrum is well approximated by the Bohr 
levels. On the other hand. Theorem O] showed that there exist true problems 



with arbitrarily small nuclear radius whose spectra depart significantly from the 
Bohr levels. In this section, we examine the true problem for higher angular 
momenta and discover that not only are the Bohr levels stable against positive 
perturbations, they are also remarkably stable against general perturbations in 
a manner we will specify below. This can be traced directly to the fact that 
C^(0, oo) is a core for the radial Hamiltonians hi, in contrast to the i = case, 
where the deficiency indices are (1, 1). 

We firstly examine the case V{r) > 0. Here, we have the analogue of Theo- 
rem pTT] : for a > we define /i^ on L'^{[a, oo), dr) by 



K=[-^.+^^j)<^' (5.1) 

where 

= I (r - a)-V e C^i[0, oo))} C L^i[a, oo), dr). (5.2) 

We have 

Theorem 5.1 For any i, let hi^tmc = he + V for some positive V G V with 
nuclear radius a > 0. Define h'f by l \5. 1\ ) for this choice of a. Then h^, /i^.tme ond 
h'(^ are bounded from below and possess infinitely many eigenvalues satisfying 

finihe) < /in(/i£,true) < fJ'niK) (5.3) 

for each n = 1, 2, . . .. 
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In the Appendix, we estimate the energy levels of /i^ as 



for n = 1, 2, . . ., where 



^"'^ (2£+l)!(2^)!(n-l)! U + W Uo/ ' ^'^ 

Thus for positive corrections of small nuclear radius, the spectrum is well approx- 
imated by the Bohr levels. 

So far we have only considered spherically symmetric potentials V. We state 
(but do not prove) an easy generalisation of Theorems |3.1| and [5.1| to the general 
case. 

Theorem 5.2 Let i^truc = H'^ + V{r) where is the idealisation with regular 
boundary conditions and V is positive and compactly supported within some radius 
a of the origin and such that r i— > 7|r|~-^ + V{r) is smooth on M^. Defining the h'^ 
for this choice of a, let 

oo 

H' = ^U*h'^U. (5.6) 

Then , Htme and H' are all hounded from below and possess infinitely many 
eigenvalues satisfying 

f^niH") < fIniHtruc) < f^n{H') (5.7) 

for each n = 1,2, 



Turning to case of general spherically symmetric potentials V & V, the sit- 
uation is as follows: the Bohr levels are preserved to arbitrary accuracy by all 
true problems of sufficiently small nuclear radius. To understand this precisely, 
consider a gedanken experiment in which one examines true problems in the 
laboratory, using a spectrometer which resolves energies down to tolerance 6E. 
Clearly, only finitely many spectral lines would be observed. In the S'-wave, one 
can find true problems of arbitrarily small nuclear radius for which the spectral 
lines depart significantly from those corresponding to the Bohr levels. In the 
sectors i > 1, however, all true problems of sufficiently small support would ap- 
parently yield the spectral lines of the Bohr levels (up to the tolerance of the 
spectrometer) with the additional possibility of extra spectral lines. The impor- 
tant point is that the structure of the Bohr levels is substantially preserved, in 
contrast to the i = case. 

These statements follow from the following result, which we prove in Section 

6. 
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Theorem 5.3 Suppose i > 1 and let A be an element of the spectrum of he. 
Given e > 0, there exists a radius a£{X, e) > such that, for any potential V G V 
supported within [0,ae{X,e)) 

dist{X,a{hi^true)) < e (5.8) 
where /if, true = he + V and dist(x, Y) = inf^gy ||s — ?/|| . 

This resuh is a direct consequence of the fact that C^(0, oo) is a core for he. 
Our statement above is an immediate coroUary of Theorem |5.3| : 



Corollary 5.4 Suppose £ >1. Given e > 0, there exists ae{e) > such that, for 
allV supported within [0,a£{e)), 

dist{\,a{he^true)) < e (5.9) 
for all X G cr{he), where he^tme = he + V. 



Proof: Given e > 0, equation (|5.9|) is non-trivial only for finitely many eigenvalues 
A of he. ae{e) may then be taken to be the minimum value of the ae{X, e) furnished 



by Theorem |5.3| over these A. □ 

An important feature of the above results is that they do not preclude the 
existence of extra energy levels, in addition to those close to Bohr levels. (This 
arises because our proof in Section 6 makes use of strong resolvent convergence.) 
Indeed, one may see that such levels can occur by considering a potential well 
within the nucleus which is gradually deepened gradually pulling bound state 
well below the Bohr spectrum. The remarkable feature of these results is that 
even during such a process, the Bohr levels remain well approximated (up to 
given accuracy) by eigenvalues of the true problem in angular momentum sectors 
i>l. 



The principal limitation of Theorem |5.3| and Corollary p.4| is that they are pure 



existence results; we have gained no information on the magnitude of ae{X,e). 
Numerical investigations, however, suggest that, for example, all true problems 
with nuclear radius a ~ 10~^ x have P-wave energy levels within 2% of each 
of the lowest five P-wave Bohr levels. 



6 Rigorous Proofs 

In this section, we provide the proofs of the theorems stated above. We begin 
by developing a sufficient condition for a sequence An of self-adjoint operators 
to converge to a self-adjoint limit A. Our condition is only slightly weaker than 
that given by Theorem VIII. 26 in [|15| (which we will also use) and the proof is 
similar. In the sequel we will often use the concepts of the graph of an operator 
and also the strong graph limit of a sequence of operators: 
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Definition 6.1 The graph T{A) of an operator A on a Hilbert space Ti is the 
set of pairs Aip) G 7^ x where ip runs through the domain of A. If {A„} 
is a sequence of operators on Ti, the strong graph hmit of the An is the set of 
pairs eHxH for which there exists a sequence 0„ satisfying 0„ G D{An), 

(pn ^ 4> and An4>n ^ V'- 

We now give our sufficient condition for strong resolvent convergence. 

Proposition 6.2 Let {An} and A be self-adjoint on Ti and let V he a core for 
A. If the graph of A restricted to V is contained in the strong graph limit of the 
An, then An A in the strong resolvent sense. 

Proof: Take (p eV. Then there exists a sequence {(pn}, (pn G D{An) with (pn ^ (p 
and AnCpn Acp. For any z G C\m, we have 

'{An + z)-^-{A + Z)-^] {A + Z)-^ = {An + Z)-^ [{A + Z)<P - {An + z)(Pn] 

+<Pn-<P (6.1) 

— > 

as n ^ oo, because the {An + z)~^ are a uniformly bounded family of operators. 
Thus we have 

{An + z)-^iP ^ {A + z)-^^P (6.2) 

for all ip G Ran((A + z) \ T>), which is dense in Ti because P is a core for A. The 
result therefore extends to all ip E Ti. because the {An + z)~^ and {A + z)~^ are 
uniformly bounded. □ 

We will also need Theorem VIII.26 of llTSi: 



Proposition 6.3 Let {An} and A be self-adjoint on Ti. An A in the strong 
resolvent sense if and only if the graph of A is the strong graph limit of the An- 

We study two classes of sequences of true problem Hamiltonians whose nuclear 
radii shrink to zero: S and S^. A particularly useful feature of such sequences is 
that, for any smooth wavefunction ip compactly supported away from the origin, 
Hnip — > (—A + 'y/r)ip as n — >• cxd, because the support of the correction to the 
Coulomb field eventually lies closer to the origin than the support of ip. In terms 
of graphs, this shows that the graph of -ffideai = —A + 7/r on C^(K^\{0}) is 
contained in the strong graph limit of any sequence in S. 

In the results below, it will be useful to have available explicit cores for 
the self-adjoint extensions of -ffideai- To this end, we define the mollified zero 
energy wavefunction of by 

x''{r) = g{r)r-'{(P{r)-Le{r)) (6.3) 

where g{r) is a smooth mollifier satisfying g{r) = 1 for r < oq and vanishing for 
r > 2aQ. Clearly ^ L^(R^,(i^r). We now define 

= C^{r\{0}) + {Xx^ I A G C} (6.4) 

and show that it is a core for H^. 
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Lemma 6.4 is a core for . 

Proof: \ is clearly symmetric and extends i^idcah so H*^^^^ extends {H^ \ 
T>^)* . Therefore the deficiency subspaces of \ T>^ are contained within those 
of -ffideai- However, a simple integration by parts argument shows that neither of 
the functions "ipQ = r~^u^ are in the domain of {H^ \ V^)*, where are defined 
in (|2.8| ). Since i/jq span the two 1-dimensional deficiency subspaces of -ffideai, it is 
clear that the deficiency indices of on are (0, 0), so that is a core for 
H^. □ 



We begin with Theorem which demonstrates the connection between the 



limit points of sequences in S and the self-adjoint extensions of -ffideai- 



Theorem |23] Let H^^e^i = -A + 7/r on C^{R^\{0}). (a) Let G S. // 



{Hn} has a self-adjoint limit H in the strong resolvent sense, then H is a self- 
adjoint extension o/ifideai- (b) Furthermore, all self-adjoint extensions o/i^ideai 
arise as the strong resolvent limits of sequences in E. 



Proof: (a) Suppose that {Hn} has a self-adjoint limit. Then by Proposition [O] , 
the strong graph limit of the Hn is equal to the graph of H. Furthermore, our 
observation above shows that the graph of H necessarily contains the graph of 
-f^ideai- Hence if is a self-adjoint extension of i^ideai- 

(b) Now suppose that is the self-adjoint extension of i^ideai with scattering 
length L. By Lemma |6.4| , we know that is a core for H^. Proposition |0 



shows us that it it is enough to construct a sequence {Hn} G S whose strong 
graph limit contains the graph of H^ restricted to T>^. In fact, as we know that 
the strong graph hmit of any sequence in S contains the graph of ifideai, it is 
enough to construct such a sequence whose strong graph limit contains the pair 
{x^ , H^X^) ■ III order to do this, we employ similar arguments to those used 



in the proof of Theorem |4.1| to construct a sequence of potentials of decreasing 
support all of which have scattering length L. Specifically, let the Xni^) be smooth 
spherically symmetric compactly supported functions on satisfying 

1. Xn{r) = X^(^) for > Cio/n 

2- Xn ^ X'^ in L'^{R^, d^r) 

3- Xn{r) is bounded away from zero for r < ao/n. 
We then define the potentials V^(r) by 

T/ _ / Xn'AXn - 7A for r < . . 

^"^"^^"1 forr>ao. ^^'^^ 

It is easy to see that the Vn & V and so {Hn} given by iJ„ = — A + 7/r + Ki(r) 
is an element of E. Furthermore, each term in the sequence has scattering length 
L. It is clear that HnXn = H^x^ fo^' ^ill n (for within radius oq, both sides 
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vanish, whilst for r > oq, Xn and are equal, and Hn and have the same 
action). Moreover, because Xn ^ i the strong graph limit of the contains 
{x^, H^X^)- By our remarks above, this suffices to prove that Hn — > if in the 
strong resolvent sense. □ 

Next, we consider the situation when we restrict from S to S+ and prove 



Theorem p.2| . (We defer the proof of Theorem |3.1| to the end of this section.) 



Theorem |3^ All sequences in converge in the strong resolvent sense to , 
the self-adjoint extension of H^^^^x with regular boundary conditions. 

Proof: Let Hn be a sequence in In exactly the same way as in the proof of 
Theorem p?T|(b), to show that Hn — >■ if", it is enough to show that {x^,H^x^) 
is contained in the strong graph limit of the Hn, where x° is the mollified zero 
energy wavefunction of H^. To prove this, let ipn be the generalised zero energy 
wavefunctions for the Hn, normalised such that ipn{f) = r~^{(l){r) — LnO{r)) for 
r > ttn, where is the radius of support and L„ the scattering length of Vn- 
Then Xn{f) = g{f)'4'n is in L^(K^,(i^r). 
Moreover, Xn ^ D{Hn) for all n and 

HY - HnXn = LnT-' {g" Q + 2g' &) ^[r - a„) (6.6) 

where dix) is the Heaviside function. Thus, because the term in parentheses 
is smooth and compactly supported, this tends to zero, since, by Theorem [4.2| , 
< L„ < 0(a„)/^(a„) ^ 0. Thus HnXn ^ if°X°- 

It therefore remains to show that Xn ~^ We know that — Xn = 
Lng{r)r~^6{r) for r > an and that g{r)r~^6{r) has finite norm. Hence 

\x\r) - Xnir)\^r^dr (6.7) 

as 77, — > oo. Furthermore, for positive potentials, we know from the proof of 
Theorem that rxn{f ) is increasing in r on (0, an) for sufficiently small a„, 
which implies that |rx„(r)| is bounded above by 0(an) — LnO{an) on this interval 
and hence that Jq" \xn - X^?r'^dr < /o""(|rx(r)P + \rXn{r)\'^)dr 0. 

We therefore have Xn ~^ which entails that H^X^) is in the strong 
graph limit of the Hn, yielding the required result. □ 

As we noted in Section 2, whilst motivating the choice of strong resolvent 
convergence, the spectrum cannot expand in the limit of strong resolvent conver- 
gence. The precise statement of this result is: 

Proposition 6.5 (Theorem VIII. 24(a) in |1^) Let {An} and A he self-adjoint 
and suppose An A in the strong resolvent sense. Then if \ ^ cr(^); there exist 
Xn € cr(y4„) such that A„ ^ A. 



This proposition allows us to prove Theorem |5T3 . 
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Theorem |5.3| Suppose i > 1 and let X be an element of the spectrum of hi 



Given e > 0, there exists a radius a£{X, e) > such that, for any potential V G V 
supported within [0,ae{X,e)) 

dist{X, cr{he^true)) < e (6.8) 
where /i£,true = hi + V and dist(x, Y) = inf^gy — y\\ . 

Proof: Suppose otherwise. Then one could construct a sequence V^(r) of poten- 
tials in V with Vn compactly supported within [0, ao/n], and hence a sequence of 
operators he^n = he + Vn{r) such that dist(A, a{hi^n)) > e for all n. But it is clear 
that the strong graph limit of this sequence contains the graph of hi restricted to 
C^(0, oo), on which it is essentially self-adjoint. By Proposition |6.2| , therefore, 
we have hi^n he. Proposition |6.5| shows that there must therefore exist a se- 
quence {A„}, with A„ G <y{he^n) and A„ — > A. But this contradicts the fact that 
dist(A, a{hi^n)) > e for all n by construction. □ 



We conclude with the proof of Theorem |5.1| , which includes Theorem ^]T| as 
a special case. This requires some more definitions. We write Q{A) for the 
quadratic form domain of an operator A and employ the usual abuse of notation 
by writing (0 | Aip) for the value of the quadratic form associated to A acting 
on 0, V'- By fin{A), we denote (if it exists) the n'th eigenvalue of a self-adjoint 
operator A which is bounded from below, counting in increasing order and with 
multiplicity. 

We shall use the min-max principle (Theorem XIII.2 in [|l^) to determine the 
eigenvalues of the operators at hand. Specifically, for A a self-adjoint operator 
bounded from below, define 

Un{A) = sup inf {ip I A^) (6.9) 

ipi,...,ip„-l Ve!7A('Pl>---.Vn-l) 

with 

Ua{^1, ...,iPm) = W\ \m=l,i^e Q{A),ij G [<^l, . . . , ipm]^} (6.10) 

where [ipi, . . . , ip^] denotes the linear span of (^i, . . . , ipm- The min-max principle 
states that, for each n, either there are at least n eigenvalues fii{A), . . . , fin{A) 
of A given by fir{A) = I'riA) for r = 1, . . . ,n below the essential spectrum of 
A, or UniA) = inf(Tess(^) in which case there are at most n — 1 eigenvalues 
below the essential spectrum given by fir{A) = i^riA) for r = 1, . . . ,n — 1 and 
PmiA) = inf(Tess(^) for all m > n. Thus if A has infinite discrete spectrum 
below the essential spectrum (which holds for the operators we consider) these 
eigenvalues are given directly by the min-max principle. In order to compare the 
eigenvalues of operators on different Hilbert spaces (as required for Theorems |3.1| 
and |5.1|) we make the following definition, which is a mild extension of a definition 



in Section XIII.2 of O . 
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Definition 6.6 Let A and B he self-adjoint and hounded helow on Hilhert spaces 
Hi, 0.2 respectively, where H2 is isometrically emhedded in Hi hy J . We write 
A<Bij and only %j JQ{B) C Q{A) and {Ji) \ AJi))n^ < {^p \ B^)^^ for all 
^ljeQ{B). 

As an immediate consequence of this definition and the min-max principle, 
we have (cf. Section XIII. 2 in |17]) 

Proposition 6.7 Let A and B he self-adjoint and hounded helow on on Hilhert 
spaces Hi,H2 respectively, where H2 is isometrically emhedded in Hi hy J . If 
A < B then z/„(y4) < z/„(i?) for each n. 

Theorems p.l| and ^]T] turn on the following lemma. Recall that Jiq denotes 
Hq, the S'-wave radial Hamiltonian with regular boundary conditions, which is 
bounded below. 

Lemma 6.8 Let /i^,true = he + V for some positive V{r) G V, supported within 
[0,a), for a > 0. For each i = 0,1,2,..., the operators /i^,true; one? h'^ are 
hounded from helow and satisfy hi < /i^.tme < h'^- 

Proof: By hypothesis on the class of potentials V from which V is taken, Qihi) = 
Q{hi^tj.^e), and since V{r) > 0, we have hi^tme bounded below and hi < hi^tme (we 
are in the special case of Definition in which Hi = H2 and is the identity). 
To prove the second inequality, let JT" be the isometry J' : L^([a, 00), dr) 1— > 
L'^{W^,dr) given by 

J ipi^r) for r > a 
I for r < a. 



(6.11) 



Q{h'f ) is the set of ip such that 

dr 

and it is clear that any such ip satisfies 

2 




dr < 00 







d 
dr 




dr < 00 



as J'lp is continuous and piecewise once continuously differentiable on 
JQ{h'^) C Q{hi) = Qihi^trne)- Clearly we have {J^p \ hi^tmeJi^) = 
all ip ^ Q{K)y so K is bounded below and /i£,true < K- n 
We now prove Theorem |0 . 



(6.12) 



(6.13) 

l+. Thus 

h'^ip) for 



Theorem |5.1| For any £, let he^tvue = hf + V for some positive V & V with 
nuclear radius a > 0. Define /i^ hy (\5.1\ ) for this choice of a. Then hi, hi^tme one? 
h'^ are hounded from helow and possess infinitely many eigenvalues satisfying 

Urihi) < /ir(Vtrue) < l^'riK) (6.14) 

for each n = 1, 2, . . .. 
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Proof: From Lemma |6.8| and Proposition |6.7| , we have Vnihi) < Vnihi^truc) < 
Pnih'i). hi has infinite discrete spectrum below its essential spectrum and in the 
Appendix, we show that the same is true for h'^, so the min-max principle gives 
l^n{h() = ^nihi) and fJ,n{h'i) = fnih'f)- A further application of the min-max 
principle shows that /i^,true also has infinite discrete spectrum below its essential 
spectrum, thus yielding the required inequahties. □ 



7 Conclusion 

We first briefiy consider the relation of the present paper to our work with Kay on 
the large scale effects of large objects in quantum mechanics [Q. More discussion 
is given in 0. The 'principle of sensitivity' states that the fine details of a small 
object in the true problem (here the specfic nuclear model chosen) can only change 
the large scale behaviour of the system if a family of idealised problems with 
differing large scale behaviour can be found. In the case in hand, we have seen in 
Section 2 that there is such a family of idealised problems in the 5'-wave, whilst 
for angular momenta ^ > 1, there is a unique idealised problem in each sector. 
Furthermore, Theorem [4.1| shows that true problems of arbitrarily nuclear radius 
can have spectra differing markedly from the Bohr levels in the S'-wave, so there 
is good agreement with the principle of sensitivity in this sector. The situation 
for £ > 1 is more subtle. We have seen (Corollary to Theorem ^.3|) that 



all true problems of sufficiently small nuclear radius preserve the Bohr levels to 
arbitrary accuracy; however, additional energy levels (termed 'rogue eigenvalues' 
in [^) may also be introduced. There is, therefore, a sense in which the spectrum 
of the idealised problem is highly stable against general perturbations in these 
sectors, in accordance with the principle of sensitivity. Our considerations have 
thus demonstrated broad agreement with, and have clarified the scope of, the 
principle of sensitivity. 

Our results in Section 5 also lead to an interesting viewpoint on the phe- 
nomenon of 'accidental degeneracy' |0, ^, ^ |T0|; i.e. the fact that the Bohr levels 
for £ > 1 are coincident with energy levels in the S'-wave. This phenomenon is 
due to the conservation of the Laplace-Runge-Lenz vector, a symmetry of the ide- 
alised Hamiltonian with regular boundary conditions which is not shared by any 
true problem. If the nucleus is modelled by a positive correction to the Coulomb 
potential of small nuclear radius. Theorems ^Tl] and |5.1| demonstrate that the de- 
generacy is approximately maintained. More generally one might expect that the 
degeneracy would be badly broken. Whilst the S'-wave energy levels can indeed 
depart significantly from the Bohr levels. Theorem ^]3| provides the surprising 
result that degeneracy is maintained to arbitrary accuracy for the eigenvalues in 
i > 1 for all true problems of sufficiently small nuclear radius, modulo the pres- 
ence of rogue eigenvalues. For angular momentum higher than zero, accidental 
degneracy is no accident! 
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Our justification of tlie regular boundary condition was premised on tlie as- 
sumptions tliat tlie nuclear structure is spherically symmetric and can be mod- 
elled by a positive correction of compact support to the Coulomb potential. Theo- 
rem |5.2| indicates how the requirement of spherical symmetry may be relaxed. We 
also note that many non-everywhere-positive potentials have scattering lengths 
L < a (for instance, if the potential is 'only slightly' negative) and that one could 
therefore use the scattering length formalism of Section 4 to estimate the energy 
levels as being close to the Bohr levels. However to do this would necessitate a 
careful study of 'believabihty'. It would be interesting to make rigorous state- 
ments about such cases and also to relax the requirement of compact support to 
allow, for example, exponentially decaying tails outside the nuclear radius. 

To summarise: by identifying the different roles of true and idealised problems 
in quantum mechanics, we have clarified the issue of which boundary conditions 
should be imposed at the Coulomb singularity of the ideahsed hydrogen atom. 
Under mild physical restrictions on the true problem, we have provided a rigorous 
justification for the regular boundary condition. Furthermore, we have developed 
a formalism for matching general true problems to (potentially) irregular bound- 
ary conditions, and have demonstrated that the entire 1-parameter family of 
boundary are required to model the full class of true problems, which (for ar- 
bitrarily small nuclear radius) can have S'-wave energy levels differing markedly 
from the Bohr levels. 
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Appendix: Spectral Properties of h'^ 

In this appendix we prove the following statements concerning the spectrum of /i^ 
on L'^{[a,oo),dr) (recall that fin{H) denotes the n'th eigenvalue of H, counting 
in increasing order and with multiplicity, where H is self-adjoint and bounded 
below) : 

1. h'g has infinite discrete spectrum below its essential spectrum [0, oo). More- 
over, /in(/if) e [fln(he),fin+lihe)). 

2. In i = 0, /in(^o) — /^n+i(^o™'"') where Lmax = 4>{o)/^{<^)i provided a < ai, 
where ai is the first non-trivial zero of 9{r) away from the origin. 
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3. For a -C Oq, 



for 



(2^+l)!(2£)!(n-l)! U + W Uo/ ' ^ 
To prove these results, we first write u^ir) for the solution to 



on such that u{r) is locally square integrable with measure dr as r ^ oo, i.e., 
u{r) oc >V|^|/2K;f+i (2'^'")- Clearly, h'^ has eigenvalue — if and only if ^^(a) = 0. 
Standard arguments give 

u^,^{a)u'^^{a) - u^^{a)u'^^{a) = {nj - nj) / u^^{r)u,,.-,{a)dr (A.4) 

J a 

(the integral converges due to the boundary conditions on 'Ufj(r)). From this it 
follows that firstly, ^^(a) has isolated zeros as a function of k, and secondly 



d ( <(r) 



^^''^^'^"^r (A.5) 



Hence u'^/uK\r=a is increasing as a function of between isolated poles, which 
correspond to eigenvalues of h'^. Now write for the solution to ( |A.3| ) with 
regular boundary conditions at r = 0. Clearly, has eigenvalues if and only if 

Similar arguments to those above show that v'i^/vn\r=a 
is decreasing but everywhere positive as a function of /t"^. We conclude that 
between every two eigenvalues of hi, there is exactly one pole of u'^/u^^^r=a^ and 
hence exactly one eigenvalue of h'^. Thus in particular, there are infinitely many 
negative eigenvalues. 

Next, we know from Lemma |6.8| and Proposition that fin{he) = Vn{he) < 
^nih'i), where i^niK) is determined by the min-max principle (see Section 6). By 
the min-max principle, I'nih'f) is either the n'th eigenvalue below the essential 
spectrum, or the infimum of the essential spectrum. Hence mi aess{h'() > 0. 
However, if i'n{h'() = inf (Tess(/i^) the min-max principle states that there are at 
most n — 1 eigenvalues below inf cress(^^)- However, we have seen above that there 
are infinitely many negative eigenvalues, so Vn{h'() is the n'th eigenvalue Unihe). 
Putting this with our earlier observation shows that lin{h'() G [^n(ht) , lin+i(ht)) . 
That aess{h'() = [0, oo) may be seen by an explicit consideration of the generalised 
eigenfunctions. 

To prove the second statement, we define tlIk) to be the position of the first 
zero (away from r = 0) of the generalised eigenfunction of at energy — k^. 
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By arguments similar to the Stiirm oscillation lemma one can see that tlIk) 
is increasing continuously in k. Moreover, ro(0) = Oi, so ro(/t) > a for all k. 
Inverting, we define L{k,) to be such that rL(K)(/t) = a. Then L{k) is decreasing in 
n where it is defined, and, by the fitting formula (|4.7| ), L(0) = Lmax = (p{ci)/9{a), 
which is strictly positive for a < ai. Moreover, because ro(/t) > a for all k, L{k) 
cannot vanish and so remains in (0,Ljnax]- It is clear that h'^ has an eigenvalue 
at — if and only if Hq^'^^ does, and also that ^inih'o) = fin+i{h^^'^^) because 
/^n+i(^o^'"'') is the eigenvalue of h^'^'^^ between fin{ho) and jJ.n+i{ho)- But we also 
know that Hm{hQ^'^^) < /iml^o""""), so the second statement is proved. 

To prove the third statement, note that it is enough to determine those values 
of K for which ^^(a) above vanishes. Now ^^(a) may be written 

u^(a) = e-'^"(2Ka)^+^^(£ + 1 + 7/2k, 2i + 2; 2Ka) (A.6) 



where '^{a,c;z) is the irregular confiuent hypergeometric function [jT8|, so it is 
enough to find the zeros of \E'(£ + 1 + 7/2k, 2i + 2; 2K,a). We look for solutions in 
the neighbourhood of the Bohr levels —7/2^ = n + i for n = 1, 2, . . .. Writing 
—7/2^ = n + i + and expanding in powers of ^, we find 

1 + 7/2K,2^ + 2;2Ka) = ^(1 - n, 2£ + 2; 2(n + £)-^a/ao) 

+K{i,n,a/ao)^ + 0{^^) (A.7) 

where 

K{i,n,S) = -^'(1 -n,2£ + 2,25(n + £)-i) 

+26} -^^{2~n,2i + 3;26{n + e)-') (A.8) 

(n + £)^ 

and \l''(a, c; z) = (9/9a)^(a, c; z). Working to lowest order in a/ao, we find 
^{i+l + ^/2K,2i + 2-2Ka) ^ (-ir'j^^^ 



/ 2 \ -(2^+1) 
-a-ir-\n-m2i)\(^-) (A.9) 



and so we find the approximate solutions given in our third statement. 
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Figure Captions: 



Figure 2.1: The function G{z). 



